Ground state properties of a Tonks-Girardeau Gas in a periodic potential 
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In this paper, we investigate the ground-state properties of a bosonic Tonks-Girardeau gas confined 
in a one-dimensional periodic potential. The single-particle reduced-density matrix is computed 
numerically for systems up to N — 41 bosons. Then we are able to study the scaling behavior of the 
occupation numbers of the most dominant orbital for both commensurate and non-commensurate 
cases, which correspond to an insulating phase and a conducting phase, respectively. We find that, in 
the commensurate case, the fractional occupation decays as /o oc 1/lniV with the particle number, 
while fo oc 1/JV ' 44 in the non-commensurate case. So there is no BEC in both cases. The study of 
zero momentum peaks shows that n(0) oc N in commensurate case, which implies that all bonons 
are localized in the insulating phase. 
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I. INTRODUCTION 

With the development of optical lattices and atom chip 
traps, quasi-one-dimensional cold atom systems have 
been realized by tightly confining the particle's motion 
in two directions to zero-point oscillation [l|, 0, 9 • Mean- 
while, using the Feshbach resonance or by tuning of the 
effective mass of particles moving in a periodic potential 
[l], IH, the inter-particle scattering length can be tuned 
to almost any desired value. These progresses have led 
to experimental realizations of the one-dimensional (ID) 
exactly solvable model that describes an interacting Bose 
gas SI Hi- 

At very low temperatures and densities, a ID Bose gas 
is expected to behave of imp enetrable particles 

known as hard-core bosons [2, fiol . In particular, two 
recent experiments successfully achieved the required so- 
called Tonks-Girardeau (TG) regime and made the TG 
gas a physical reality [H,12|. Physically, a TG gas is de- 
fined to be a ID strongly correlated quantum gas con- 
sisting of bosons with hardcore interaction. This model 
of a ID gas was first proposed by Tonks in 1936 [Hj]. As 
a milestone development to the model, Girardeau solved 
the model exactly by the famous Bose-Fermi mapping 
theorem [j| H, [l3j], where the TG gas can be mapped 
to a spinless free fermion gas. Recently, it was found 
that the above case is just a special case of a general 
mapping theorem between bosons and fermions in one 
dimension, where the particles can interact with finite 
strength [ME!. 

Since the TG gas is both theoretical exactly solvable 
and experimental accessible, there are great research in- 
terest recently in the T G g as in different trapping po- 

tgntiais EMP M 0, HI, M H M M, m M 

1281 129 . [30L l3ll . l32j . A review on relevant studies in this 
interesting subject can be found in a recent article by 
Yukalov and Girardeau [l3[. In particular, Lin and Wu 
investigated the grou nd-state properties of a TG in pe- 
riodic potentials [2g|. They used the Monte Carlo inte- 
gration technique to compute the single-particle reduced 



density matrix (SPRDM) for small systems up to N = 7 
bosons. Their analysis of ground states shows that when 
the number of bosons N is commensurate with the num- 
ber of wells M in the periodic potential, the boson system 
is a Mott insulator whose energy gap, however, is given 
by the single-particle band gap of the periodic potential; 
when N is not commensurate with M, the system is a 
metal (not a superfluid). The purpose of this work is to 
compute the SPRDM for large systems containing more 
than 40 particles, and then study scaling relations, quan- 
titative predictions for scaling exponents of the ground- 
state occupation numbers, and zero-momentum peak in 
different filling conditions. Therefore we put particular 
emphasis on examining the ground-state fractional occu- 
pation in these two phases of the system. 

The remaining of this paper is organized as follows. 
In Sec. II, we introduce the model Hamiltonian of the 
system and describe the single-particle eigenstates and 
eigenvalues which will be used. In Sec. Ill, we review 
the Bose-Fermi mapping theorem and then construct the 
exact many-body ground state wave function. In Sec. IV, 
we devote ourselves to study the many-body properties 
of the TG gas. Summary and conclusion are given in 
Sec. V. 



II. MODEL HAMILTONIAN AND 
SINGLE-PARTICLE EIGENSTATES 



Model Hamiltonian 



We consider a gas of N hardcore bosons trapped in a 
tight atomic waveguide. The waveguide restricts strongly 
the dynamics of the gas in the transversal directions, such 
that in the low temperature limit we can define our model 
in the longitudinal direction only. In this direction we 
consider a periodic potential such that the many-particle 
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Hamiltonian at low density can be written as 



N 



H = J2 



2m dx 



(1) 



where m is the mass of a single boson, V(x) = V(x + 
d) the periodic potential with d being the period. In 
this work, we use the Kronig-Penney (KP) potential as 
a candidate for the periodic potential, which takes the 
form 



M 



V(x)=jJ2 6 ( x -^ 

3=1 



(2) 



where 7 is the strength of the S potential and M is the 
total number of wells. 



III. BOSE-FERMI MAPPING THEOREM AND 
MANY-BODY WAVE FUNCTIONS OF THE TG 
GAS 

For a system of N identical hardcore bosons in the ID 
external potential V(x), the bosons interact with each 
other by impenetrable pointlike interactions, which can 
be more conveniently treated as a boundary condition for 
the many-body wave function \& b{x±, x%, . . . , Xn, t): 

\1/ b (xi , X2 , ■ ■ ■ , xn, t) = if Xi = Xj , 1 < i < j < N. (9) 

Then the hardcore boson gas can be considered as a 
free boson gas governed by the following free Schrodinger 
equation, 



d 



N 



n 2 d 2 



(10) 



B. Eigenstates and eigenvalues of periodic 
potential 

We impose the usual periodic boundary conditions and 
use d as the unit of distance and h 2 /(2m) as the unit of 
energy. Then we arrive at the following dimensionless 
single-particle Hamiltonian, 



3=1 



(3) 



The Hamiltonian can be exactly solved. One can find the 
solution in any standard textbook of solid state physics. 
We review the result again for later computation. The 
Bloch wave functions take the form 

il>a{x) — C a [sm(k a x) + e~ tKa sin[/c a (l — x)]], < x < 1, 
ip a (x) = e iK "Wip a {x -[»]), 1 < x < M. (4) 

The eigenenergy is 

E a = k 2 a (5) 
where k a satisfies the transcendental equation 
sin(fc Q ) 



cos(fc Q ) 



T- 



cos(K a ), 



(0) 



2ira M - 1 . , 

K a = -jj-, a = 0,±l±2,---±— ^— (7) 

and the normalization constant for the Bloch wave func- 
tion is given by 



CVy — 



M(k a — i sin(2fc Q ) + cos K a (sin k a — k a cos k a ) 



(8) 



where the wave function ^ b satisfies the hardcore bound- 
ary condition of Eq. (9). Based on the observation that 
the hardcore boundary condition of Eq. (9) is automat- 
ically satisfied by a wave function of fermions due to its 
antisymmetry, Girardeau gave the exact many-body 
wave function of the hardcore boson system via the fa- 
mous Bose-Fermi mapping, which relates the wave func- 
tion of hardcore bosons to that of noninteracting spinlcss 
fermions in the same trapping potential: 



Vb(xi,X2, ■ ■ -,x N ) = A^ F (xi,x 2 , ■ ■ ■ ,x N ,t). (11) 



with 
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Yl sgn(xi 

l<i<j<N 



(12) 



where sgn is sign function and A is unit antisymmetric 
function which ensures that ^b has proper symmetry 
under the exchange of two bosons. The free fermionic 
wave function can be compactly written in form of the 
Slater determinant 



1 



=Det 



N 



-A^m(Xj,t)]. 



(13) 



where ip m ,m = 1, ...,N is single-particle eigenstate. 
The eigenstates are governed by a set of uncoupled single- 
particle Schrodinger equations, 



ih^-ip m (x) 
at 



2m dx 2 



tp m (x). 



(14) 



Eqs. (11), (13) and (14) describe how to construct the 
exact many-body wave functions of a TG gas in any ex- 
ternal potential V(x). Then the wave function of the TG 
gas is 



A(xi,x 2 , ■ ■ - ,x N )- 



^Det^ , 

/]vT J ' 



,lbPm{Xj,t)]. (15) 
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where ip m (x) is the single-particle eigenstate in the po- 
tential V{x). 

For the periodic potential that we are interested in, 
the single-particle eigen function is given by Eq. (4), then 
we can construct the exact many-body wave function of 
the TG gas in the KP potential according to the above 
procedures. The ground-state many-body properties of 
the TG gas in the KP potential can be studied from this 
exact wave function. 
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FIG. 1: Single-particle reduced density matrix, p(x, x'), of the 
TG gas in periodic potential for different particle numbers, 
(a) AT = 1, (b)N = 15,(c)iV = 25, (D)JV = 35 when number of 
wells M = 35. 



IV. MANY-BODY PROPERTIES OF TG GAS IN 
PERIODIC POTENTIALS 

In this section we numerically evaluate the ground- 
state properties of the ID TG gas in the periodic po- 
tential using the exact many-body wave function of the 
previous section. 



A. Single-particle reduced density matrix 

The many-body wave function \&s fully describes the 
state of the system. However, its form does not transpar- 
ently yield physical information related to many impor- 
tant observable such as occupation numbers of natural 
orbital and the momentum distribution. The expecta- 
tion value of one-body observable are readily obtained 



from the SPRDM, 

p(x,x') = N J dx 2 ■ ■ ■ x N ^* B (x,x 2 , ■ ■ ■ ,xn) 

x$ b (x',X2,...,x n ), (16) 

which satisfies the normalization condition 



p(x, x)dx = N 



(17) 



Although the exact many-body wave function of the 
TG gas can be written in a compact form, the calcula- 
tion of the SPRDM is a very difficult task as it is very 
time consuming to calculate multidimensional i nteg rals 
in SPRDM for a large system [fl [H, M, HH HI M, HI 
H! OH li. However, in the TG limit the SPRDM can 
be written in a matrix product form in terms of single- 
particle eigenstates |29j 



p(x,x') 



N 

E 



^*(x)A lJ (x,x')ip J {x'). 



(18) 



Here the N x N matrix A(x, x') = {Aij (x, x')} takes the 
form, 



(19) 



A = (p-y Det P 

where the entries of the matrix P are Pij(x, x') = 5ij — 

2 if)*(t)ipj(t)dt, and we can assume x < x' without 
loss of generality. This formalism of the SPRDM enables 
us to calculate considerable large systems of the TG gas. 

The SPRDM expresses self-correlation and one can 
view p(x,x') as the probability that, having detected 
the particle at position x, a second measurement, imme- 
diately following the first, will find the particle at the 
point x'. Classically, p(x,x') = S(x — x'), so the off 
diagonal elements of SPRDM comes from purely quan- 
tum correlations of the particles. Fig. Q] displays con- 
tour plots of SPRDM, p(x, x'), for the number of particle 
N = 1,15,25 and 35 respectively, where the number of 
wells of the periodic potential M = 35. We clearly see a 
characteristic pattern for each value of iV: The SPRDM 
are largest close to the diagonal elements which stands for 
the position density distributions. The diagonal elements 
of SPRDM shows oscillations due to the barrier of the 
periodic potential which tends to repel the particles and 
push the particles stay inside the well. The off-diagonal 
elements of the SPRDM relate to off-diagonal long range 
order (ODLRO) [36j and we can see that the off-diagonal 
elements are decreasing in contrast to the diagonal as 
the number of particles N increases. It means that the 
repulsive interaction tends to destroy the ODLRO, and 
ODLRO vanishes for a system of hard core bosons in a 
ID periodic potentials in the thermodynamic limit. As 
N = M, the SPRDM is almost localized, and the system 
becomes an insulator. 
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FIG. 2: Fractional occupations, fo = Xo/N, as a function of 
particle number N at commensurate case, N/M = 1, when 
periodic potential strength 7 = 2. 
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FIG. 3: Occupation numbers of lowest orbital Ao, as a func- 
tion of particle number, l/ln(iV), at commensurate case, 
N/M — 1, when periodic potential strength 7 = 2. 



B. Occupation numbers and natural orbital 

In a macroscopic system, the presence or absence of 
BEC is determined by the behavior of p(x,x') 
ODLRO is present if the largest eigenvalue of p(x,x') is 
macroscopically non- vanishing. In this case the system 
exhibits BEC and the corresponding eigenfunction, the 
condensate orbital, plays the role of an order parame- 
ter [35|, [36|. The fraction of particles that are in the 
lowest orbital is related to the largest eigenvalue of the 
SPRDM by f = \ /N. Therefore, in analogy to the 
macroscopic occupation of a single eigenstate in a Bose- 
Einstein condensate, this orbital is sometimes referred 



to as the "BEC state and the quantity hence acts as a 
measure of the coherence in the system. 

The occupation numbers of the natural orbital are de- 
fined as eigenvalues of the SPRDM, 



dx'p(x,x')(f>i(x') = \i<f>i(x),i = 0,1, 



(20) 



where <fii(x) are the so called natural orbitals, which are 
obtained from the eigenfunctions of the SPRDM and it 
represents an effective single-particle state. The corre- 
sponding eigenvalue A^ is occupation number of the i-th 
natural orbital and satisfies ^ A^ = N. The SPRDM is 
diagonal in the natural orbital basis. The corresponding 
eigenvalue of each orbital gives the population probabil- 
ity of that orbital. The natural orbitals can be ordered 
according to their eigenvalues, we assume Ao > Ai > ... 
without loss of generality. 
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FIG. 4: Fractional occupations, fo — \o/N, as a function of 
particle number N, at non-commensurate case, N/M = 1/3, 
when periodic potential strength 7 = 2. 

The importance of natural orbitals is that, for a boson 
system, one can infer whether the system is a BEC state 
or not from scaling behavior of the occupation number 
of the lowest natural orbital Ao [H, [36[. The system 
will show a BEC behavior if the occupation number of 
the lowest natural orbital is proportional to the particle 
number of the system, which means a macroscopic num- 
ber of particles will condensed to the lowest orbital. Then 
the corresponding natural orbital will be the eigenfunc- 
tion of the BEC. For the TG gas in periodic potentials, 
there are two different phases in the ground state [28| : 
one is a Mott-insulator for the commensurate case where 
N/M is an integer. In this phase, one or more Bloch 
bands have been fully occupied. The other is a boson 
conductor phase for non-commensurate case where N/M 
is a fractional number. In this phase, the Bloch bands 
are partially occupied. What we are interested here is 
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FIG. 5: Occupation numbers of lowest orbital in(Ao), as a 
function of particle number, ln(iV), at non-commensurate 
case, N/M = 1/3, when periodic potential strength 7 = 2. 

the BEC probability and the scaling behavior of the oc- 
cupation number of lowest orbital in different phases. Di- 
agonalizing the SPRDM numerically, we can obtain the 
occupation number of lowest orbital for different particle 
numbers up to 41. 

We have shown the fractional occupations of the low- 
est orbital, fo = \q/N, with particle numbers N, at com- 
mensurate case (the lowest Bloch band has been fully oc- 
cupied) and 7 — 2 in Fig. [5] From the figure, we can see 
that the fractional occupations of lowest orbital decrease 
with the particle number increases and it will vanish in 
the thermodynamic limit. By doing finite size scaling in 
Fig. [31 we found that the fractional occupation of lowest 
orbital decreases with the system sizes as fo oc l/ln(7V) 
in the Mott-insulating phase and it vanishes in the ther- 
modynamic limit. So we can infer that there is also no 
BEC for the commensurate case. 

The fractional occupations of the lowest orbital, fo = 
Xo/N, with particle numbers N, in non-commensurate 
case for N/M = 1/3 and 7 = 2 is shown in Fig. [4] We 
can see that the fractional occupations of lowest orbital 
decrease with the particle number increases and will van- 
ish in the thermodynamic limit. Doing finite size scaling 
of occupation number of lowest orbital in Fig. [FJ from 
which we can find that the occupation numbers of lowest 
orbital Ao shows power law dependence on the particle 
number ./V of the system in non-commensurate case, i.e 
Ao ~ TV 56 . So there is no BEC behavior also in the 
non-commensurate case. 



C. Momentum density distribution 

Although the position density profiles and energy are 
exactly the same between the hardcore boson gas and 



FIG. 6: Momentum distributions for several particle number 
at commensurate case, N/M = 1, when periodic potential 
strength 7 = 2. 
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FIG. 7: n(0), as a function of particle number N at com- 
mensurate case, N/M = 1, when periodic potential strength 
7 = 2. 

their corresponding spinless fermion conterpart due to 
Bose-Fermi mapping theorem, the momentum distri- 
butions differ considerably from each other [H, |32| . 
The momentum distribution can be obtained from the 
SPRDM as 

n(Jb) = 27T- 1 J dxdx'p{x, x')e- lk{x - x,) . (21) 

The momentum distributions of the TG gas in the com- 
mensurate case are shown in Fig. [6] for various particle 
numbers. We can see that the momentum distributions 
of the TG gas has a bosonic structure, which has peaks 
at k = and the profiles are narrow. Although a TG gas 
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FIG. 8: Momentum distributions for several particle num- 
ber at non-commensurate case, N/M = 1/3, when periodic 
potential strength 7 = 2. 



has the same energy spectrum and position profiles with 
its fermionic conterpart, their momentum distributions 
are very different. In addition, the profiles of momentum 
distributions is nearly unchanged with the particle num- 
ber increases at commensurate case. This is because in 
the Mott-insulating phase the particles are localized in 
the real space. The oscillations in the momentum dis- 
tribution comes from the finite size effect, it will vanish 
in thermodynamic limit. What is more interesting here 
is that we find the momentum distributions have same 
zero momentum peaks after rescale. Then we plot n(0) 
as a function of particle number in Fig. [7] It clearly 
shows that n(0) oc N in the commensurate case. Simi- 
lar behavior has been found for the TG gas in the har- 
monic trap [24J. Thus, the system of hard-core bosons 
mimics the macroscopic occupation of a momentum zero 
state and, in this aspect, resembles a uniform and non- 
interacting Bose system. The momentum distributions 
of non-commensurate case (fractional filling) are shown 
in Fig. [5] for various particle numbers. We can see that, 
in the fractional filling case, the profiles of momentum 



distribution become narrower with the particle number 
increases, which is different from that in the commensu- 
rate case as we discussed before. This is because in this 
case, the system is a boson conductor, so the spatial dis- 
tribution become more uniform with the particle number 
increases due to many body repulsions. The small oscil- 
lations in the momentum profiles come from finite size 
effects and it will vanish in the thermodynamic limit. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the ground-state properties of 
a ID TG gas in a periodic potential. Based on the ex- 
act many-body wave function, obtained from Bose-Fermi 
mapping theorem, we have calculated the SPRDM for 
considerable large system by employing the technique 
that the SPRDM can be expressed in a matrix product 
form in the TG gas limit. We calculated the occupa- 
tion number of the lowest orbital by diagonalizing the 
SPRDM numerically. We found that fractional occupa- 
tion of the lowest orbital decreases with particle num- 
ber increases as /o oc l/\n(N) in the commensurate case 
while as fo oc 1/N in the non-commensurate case. So 
we conclude that there is no BEC for the TG gas in pe- 
riodic potential in both the commensurate and the non- 
commensurate cases. We suspect that the scaling expo- 
nent of fractional occupations in the non-commensurate 
case may approach to 0.5 in thermodynamic limit. How- 
ever, lack of data on bigger lattices prevented us to per- 
form extensive finite size scaling analysis. We may seek 
analytic approaches to derive the observed scaling behav- 
ior of lowest orbital with particle number in the future. 
The momentum distributions of the TG gas show bosonic 
behavior and zero momentum peaks increase with the 
particle number as n(0) oc N at the commensurate case. 
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